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Aims. We derive the value of the dark matter density at the Sun's location (p Q ) without globally mass-modeling the Galaxy. 
Methods. The proposed method relies on the local equation of centrifugal equilibrium and is independent of i) the shape of the dark 
matter density profile, ii) knowledge of the rotation curve from the galaxy center out to the virial radius, and iii) the uncertainties and 
the non-uniqueness of the bulge/disk/dark halo mass decomposition. 

Results. The result can be obtained in analytic form and it explicitly includes the dependence on the relevant observational 
quantities and takes their uncertainties into account. By adopting the reference, state-of-the-art values for these, we find p B = 
0.43(1 1)(10) GeV/cm 3 , where the quoted uncertainties are respectively due to the uncertainty in the slope of the circular-velocity 
at the Sun location and the ratio between this radius and the length scale of the stellar exponential thin disk. 

Conclusions. We obtained a reliable estimate of p Q > that, in addition has the merit of being ready to take into account any future 
change/improvement in the measures of the observational quantities it depends on. 
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Galaxy rotation curves (e.g. [Rubin et al.|| 1980| |Bosma et al. 
1981) have unveiled the presence of a dark "mass component" in 



spirals. They are pillars of the paradigm of massive dark halos, 
composed of a still undetected kind of matter surrounding the lu- 
minous part of galaxies. The kinematics of spirals shows univer- 
sal systematics (Persic, Salucci 1996 1 Salucci et al. 2007 1, which 
seems to be at variance with the predictions emerging from sim- 
ulations performed in the A cold dark matter (ACDM) scenario, 
|1996[ l, the currently preferred cosmolog- 



<^> (e.g. |Navarro et al. 



ical paradigm of galaxy formation (e.g. Gentile et al. 2004) 



Individual and coadded rotation curves (RCs) of spiral galaxies 
are also crucial to investigate frameworks alternative to the stan- 
dard paradigm of cold collisionless DM in Newtonian gravity 
(e.g. |Sanders, McGaugh|2002[|Berezhiani et al.|2009] ). 

At the same time, dedicated searches of DM particle can- 
didates have seen an important boost in recent years with rel- 
evant and costly experiments being planned and executed. The 
so-called direct-detection experiments look for the scattering of 
DM particles off the nuclei inside the detectors (e.g. CDMS, 
XENON 10, DAMA/LIBRA) by typically measuring the de- 
posited energy or its annual modulation. Clearly in all these ex- 
periments the signal is proportional to the DM density in the 
Sun's region, p Q . On the other hand indirect-detection experi- 
ments (in particular Super-Kamionkande, AMANDA, IceCube 
and ANTARES) search for the secondary particles (neutrinos 
in these cases) produced by DM annihilations at the center of 
the Sun or Earth, where it is expected that DM accumulates af- 
ter losing energy via scattering, possibly reaching a thermalized 
state. The expected signal in this case depends on the DM den- 
sity inside these objects, which in turn is driven, via the capture 
mechanism, by the same halo DM density in the Sun region, p G . 



Therefore, in both these kinds of direct and indirect searches, 
an estimate of the the local density p is very important for a 
precise estimate of the signal or at least reliable bounds on the 
DM cross-section vs mass to be compared with limits from other 
searches. 

What is then the value of p ? A value of 

,3 



p =0.3GeV/cm J 



(1) 



is routinely quoted in hundreds of papers, but how does this num- 
ber come out? In which works do we find the details of its mea- 
sure? It is worth observing that in most of the cases in the liter- 



ature, the above value is given with no reference (e.g. Donato et 
al.|20 09 1 Savage et al.|2009 1. Sometimes, the reference goes to a 
couple of seminal papers. Among them, the Particle Data Group 
Review ( |PDG|[2008] ) indicates the above value "within a factor 
of two or so" and justifies it as coming from "recent estimates 
based on a detailed model of our Galaxy". However, the works 
cited are neither recent nor detailed and sometimes not even an 
independent estimation of p . 

The only exception is the work by Caldwell, Ostri ker| ( |1981[ ) 
that devised what can be considered as the standard method (CO 
hereafter) to determine the value of p Q from observations (see 
below). Their resulting value, 0.23^ j2GeV/cm 3 , arises however 
from very outdated kinematical observations and from a cored 
(rather than a cuspy) halo distribution, so it is not a great support 
for equation ([TJ. Similar conclusions can be drawn by looking 
at other influential reviews: the papers they cite to back up the 
value ([T]) either do not estimate this quantity or use very outdated 
observations. 

In general, it is quite simple to infer the distribution of 
dark matter in spiral galaxies. Spiral's kinematics, in fact, reli- 
ably traces the underlying gravitational potential (Persic, Salucci 
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1996 ; Salucci et al. 2007). Then, from coadded and/or individual 
RCs, we can build suitable mass models that include stellar and 
gaseous disks along with a spherical bulge and a dark halo. More 
in detail, by carefully analyzing (high quality) circular velocity 
curves, with the help of relevant photometric and HI data, one 
can derive the halo density at any desired radius. The accuracy 
of the "measurements" is excellent and the results are at the core 



of the present debates on Galaxy formation (e.g. Gentile et al. 
[20041 |2005l |DeBlok|2009) . 

To measure p Q , instead, is far from simple, because the 
MW kinematics, unlike that of external galaxies, does not trace 
the gravitational potential straightforwardly. We do not directly 
measure the circular velocity of stars and gas but rather, at 
our best, the terminal velocity Vj of the rotating HI disk, and 



this only inside the solar circle (e.g. McClure-Griffiths, Dichey 
2007) 1. This velocity is related to the circular velocity V(r), for 
r < R Q , by means of V(r) = VrM+V© r/R B , where R Q 8 kpc is 
the distance of the Sun from the Galaxy center and V the value 
of the circular velocity at the Sun's position. Both quantities 
are known within an uncertainty of 5% - 10% (e.g. |McMillan,| 
Binney 2009) 1, which triggers a similar uncertainty in the derived 



magnitude and slope of the circular velocity. 

As a result, and also considering other kinematical observa- 
tions (see Sofue 2009), the circular velocity of the MW from 2 
to 8 kpc can be only derived within non negligible uncertainties: 



V(r) = (215 ± 30) km/s 
d\ogV(r)/dlogr = a(r) = 0.0 ± 0.15 . 



(2) 
(3) 



Note that the range of circular velocities in equation Q is cre- 
ated by a mix of a) observational errors, b) uncertainties in the 
values of R and V Q , and c) actual radial variations of V. The first 
two trigger also part of the range of the velocity slope ([3]). We 
stress that data show that the radial variations of a(r) are small 
and likely caused by the uncertainties just discussed: in either 
case we have da(r)/dr ± 0.03 /kpc 0. In other words, in 
this region the RC can be approximated by a straight line, whose 
slope is known however only within a degree of uncertainty. 

The outer (out to 60 kpc) MW "effective" circular veloc- 
ity V(r) = (GM(r)/r) l t 2 is, instead, much more uncertain and 
depends on the assumptions made on dynamical and structural 
properties of its estimators. It appears to decline with radius, 
with quite an uncertain slope (Battaglia 20 031 P^ e e * al.|2 008 i 
|Brown et al.|2009) 

dlogV(r)/d\ogr = -0.20^ (R Q < r < 60 kpc) . (4) 

These uncertainties, combined with the intrinsic "flatness" of 
the RC in the region specified above (that complicates the mass 
modeling even in the case of a high-quality RC (Tonini, Salucci 
|2004| l), make it very difficult to obtain a reliable bulge/disk/halo 
mass model and consequently an accurate estimate of p . 

To overcome these serious difficulties Caldwell and Ostriker 
( |1981| l developed a method in which other observational data, 
linked in various ways to the gravitational potential, help with 
the mass modeling. These include the l.o.s. dispersion veloc- 
ities of bright tracers at known distances from the Sun (e.g. 
OB stars) and the total Galaxy mass. This extra information 
allows a determination of p Q , though it turns out to be uncer- 
tain within a factor 2 (Caldwell, Ostriker 1981), or somewhat 
less when more constraints from the z motions of disk stars are 



added (e.g. Oiling, Merrifield 2001 Weber, de Boer 2009 Sofue 
2009)[J 



1 To overcome the observational limits of equation {3} by resorting to 
additional tracers of the MW gravitational potential also requires taking 



By averaging the results from different determinations ob- 
tained so far, one finds p = (0.3 + 0.2) GeV/cm 3 , where the un- 
certainty in the result is triggered in an unspecified way from ob- 
servational and fitting uncertainties, as from biases arising from 
the various assumptions taken by the method, not the least the 
assumed DM density profile. 

Recently, by applying the CO method with a refined sta- 
tistical analysis to a large set of observational data, |Catena. 



Ullio ( 2009| l claimed a measure with a very small uncertainty: 
p Q = (0.389 + 0.025) GeV/cm 3 . While this result woul d be no- 
ticeable, it has not been confirmed by a subsequent work (Weber, 
de Boer 2009 ) and it seems unlikely, in view of equations (|2Jl and 
Q, reflecting the state of art of our (lack of) knowledge. 

The aim of this work is to derive p by following a more 
direct route than mass modeling the whole Galaxy and dynami- 
cally modeling a number of galactic components and a series of 
galactic potential tracers. This will be done by means of a specif- 
ically devised method and by using some recent results obtained 
for external galaxies ( Salucci et al. 2007) . The idea is to resort to 
the equation of centrifugal equilibrium, holding in spiral galax- 
ies (see|FallEfttaEOTll980]for details) 



V 2 /r : 



a H + a D + a B . 



(5) 



where a H , a D , and a B are the radial accelerations generated 
by the halo, stellar disk, and bulge mass distributions. Taking 
first the (quite good) approximation of spherical DM halo, we 
have an <* r~ 2 L Ph(R) R 2 dR. A similar relation holds for the 
bulge. Therefore, by differentiating equation |5]), we obtain the 
DM density at any radius in terms of the local angular velocity 
cj(r) = V/r, the RC slope a(r), the disk-to-dynamical mass ratio 
/3(r) (see later), and the bulge mass density: 



p H (r) = aj(r) 2 [F rot (a(r)) - F D (J3(r)] - p B (r). 



(6) 



with F tot and Fp known functions. 

In spirals, equation (|6ji is not useful for determining the DM 
density at any radius because 1) it virtually collapses for r < 
R D where F tot F D , and the bulge density can also become 
dominating, p B » Ph\ 2) the radial variations of a(r) have non- 
negligible observational uncertainty that further complicates the 
effect discussed in the previous point; 3) the quantity co is known 
with less accuracy than V, the observational quantity entering the 
traditional mass modeling. 

Instead, in estimating p Q , i.e. the density of the MW DM halo 
at a specific radius (the Sun position), the above drawbacks dis- 
appear: 1) since R G > 3Rd we have F to ,(R Q ) » Fd(Rq), equa- 
tion (|6]) does not collapse and, as a bonus, the most uncertain 
term of the r.h.s. of (|6]l is also the smaller one; 2) a> Q is very 
precisely measured; 3) da/dr\/{ 0; and 4) at the Sun's po- 
sition the bulge density p B (R ) is totally negligible, < p#/50 
(e.g. Sofue et al.|2009 1. Thus, this method is very powerful for 
determining the value of the DM density at R Q . The DM density 
at any radii is obviously left to the standard mass modeling. 

The method is obviously simpler for a spherically symmet- 
ric DM halo, and can be further simplified by considering an in- 
finitesimally thin disk for the distribution of stars in the Galaxy. 
However, below we also include the effects of a possible halo 
oblateness and disk thickness. Here, we anticipate that these ef- 
fects, constrained by observations, are rather weak, of the order 
of a few percent, and are therefore irrelevant for this work. As 



extra assumptions, and the complex procedure can trigger some mild in- 
ternal inconsistency in the resulting MW mass model (see Appendixjcj. 
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a result, we obtain a reliable and model-independent determi- 
nation of the local DM halo density, as well as of its intrinsic 
uncertainty. 

In the next section we describe the method in detail and de- 
rive p Q and its uncertainty as a function of the relevant observ- 
ables. In the last section we discuss the results and draw the con- 
clusions. In the Appendices, we explicitly describe the effect of 
the halo oblateness and disk thickness, and we comment on the 
inherent problems in the traditional determination of p Q . 

2. A model-independent method 

We model the Galaxy as composed by a stellar exponential 
thin disk (Freeman |1970| l, plus an unspecified spherical DM 
halo with density profile p//(r). For the present work we can 
neglect the HI disk because its surface density, between 2 kpc 
and R , is 100 to 5 times smaller than the stellar surface den- 
sity (Nakanishi, Sofue 2003). Similarly, we neglect the stellar 
bulge because, as mentioned above, its spatial density at R Q is 



virtually zero (e.g. Sofu e et al.|2009) >. The standard method and 
its variants cannot take these very simplifying assumptions be- 
cause the global modeling of the Galaxy involves these mass 
components in a crucial way. 

As discussed in the introduction we can rewrite the equa- 
tion of centrifugal equilibrium by subtracting the disk compo- 
nent from the total acceleration. From its radial derivative we 
then find 



X„ 



4-nGr 2 dr 



r 2 I — — - a D (r) 



(7) 



where X q is a factor correcting the spherical Gauss law used 
above in case of oblateness q of the DM halo. 

We describe the factor X q in appendix [A] Since the obser- 
vational evidence for an oblateness of the Galactic DM halo is 
still quite uncertain and at the same time the mean value of q, 
as measured in a very large sample of external spirals, is very 
near to 1 (see e.g. |Q'Brien et al.| ( |2010| l for a recent review of 
this issue), we consider q = 0.95 + 0.05 here as the reference 
value and compute X q accordingly (see Appendix |A|). We thus 
find that from present observations X q boils down to a correc- 
tion of less than 5% to p , X q 1.00-1.05. Let us stress that 



by equation ( A. 1 1 one may also take into account a value of q, 
emerging from future improved observations, outside the range 
considered here. 

The disk component can be reliably modeled as a Freeman 



stellar exponential thin disk of length scale ([Picaud, Robin 2004 



M. Juric et aL||2008l |Robin et al.||2008| lReyTe||2009| l R D 



(2.5 + 0.2) kpc. The stellar surface density is then: S(r) = 
(M D /27rRl)e- r/RD . Also, the disk can be considered infinitesi- 
mally thin. In fact, its thickness zo is small, zo ~ 250 pc ( | Juric | 
|et al.||2008] l and moreover zq <k Rd < Rq, so that its effect on 
the derivative of the acceleration, and in turn on our measure, is 
very limited. For the sake of completeness, we compute it ex- 
plicitly in Appendix [B] and show that it implies a reduction of 
less than 5% of the a D term as computed for an infinitesimally 
thin disk. We can thus write a D (r) = ^^(I K - hK x )X Zo , 
where /„ and K n are the modified Bessel functions computed at 
r/2Ro, and X Zo 0.95 accounts for the nonzero disk thickness 
(See Appendix IBJ. 

Since only the first derivative of the circular velocity V(r) en- 
ters in (|7]l and in any case this function in the solar neighborhood 
is almost linear, we can write 



where a Q = a(R Q ) is the velocity slope at the Sun's radius. Then 
equation (j7]i becomes 



PH(r) 



AnG 



VHr) 



(l + 2or )- 



GMp 



-H(r/R D )X :o 



(9) 



with 2H(r/R D ) = (3I K - I X K X ) + (r/R D )(IiK - I Ki). 
Equation (|9]l holds at any radius outside the bulge region and 
measures Ph(Rq) = Po by subtracting the "effective" density of 
the stellar component from the one of all the gravitating matter. 

The disk mass can be parametrized (Persic, Sal ucci|1990[ ) by 
M D =pM.l G- l V 2 R Q , withjS= V 2 /V 2 \ Re , i .e .the fraction of the 
disc contribution to the circular velocity at the Sun. 

Finally, by exploiting the fact that the quantity V/R\r = a> = 
(30.3 + 0.3) km/s/kpc is measured with very high accuracy and 



much better than V and R Q separately (McMillan, Binney 2009 
Reid|20 09), and after defining r Q o = R Q /Ro, we obtain 



p Q = 1.2x10 



-27 



cm 3 \km/skpc 



^[(l+2a )-l.lj8/(r 0D )X a 
(10) 

where f(r @D ) = r^ D H(r QD ). 

It is now possible to observe the advantages of the proposed 
method: a) it does not require assuming a particular DM halo 
density profile, or the dynamical status of some distant tracers 
of the gravitational field; b) it is independent the (poorly known) 
values of V and of the RC slope at different radii; c) it does 
not depend on the structural properties of the bulge, which in the 
mass modeling creates a degeneration with the stellar disk and 
DM halo, d) it only mildly depends on the ratio r Q o, as well as 
on the disk mass parameter /3; finally, e) the method depends on 
the RC slope at the Sun a , although in a specified way. In all 
points a) - e) the method brings an evident improvement over 
the CO one. 

To proceed further we discuss the parameters appearing in 
equation ( 10 1. Our determination does not depend on the value 



of R Q and Rd separately, but only on their ratio r e o- For this we 
adopt the reference value and uncertainty r e o = 3.4 + 0.5^] as 
suggested by the values of R D given above and by the average of 
values of R Q obtained in recent work R Q = 8.2 ± 0.5 kpc (jGhez 
eT"aT[ ( [20081 ); |Gillessen et aT] ( [20081 ); |Bovy et al.| pOO^) ). This 
leads to f(r QD ) ^ 0.42 + 0.20, whose uncertainty propagates 
only mildly into the determination of p Q , because the second 



term of the r.h.s. of equation ( 10 1 is only one third of the value of 
the first. (And it can only reach one half by stretching all other 
uncertainties.) 

Present data constrain the slope of the circular velocity at 
the Sun to a central value of a Q — and within a fairly narrow 
range -0.075 < a Q < 0.075. The uncertainty of a Q is the main 
source of the uncertainty of the present determination of p , and 
let us recall that also values outside our adopted range may be 
used in the present analytic determination, for instance, a value 
belonging to a wider range allowed for a Q claimed by (Oiling, 
|Memfield|200l) . 

In equation ( 10 », /J is the only quantity that is not observed 



and therefore intrinsically uncertain. We can, however, constrain 
it by computing the maximum value /3 M for which the disk 
contribution at 2.2 Rq (where it has its maximum) totally ac- 
counts for the circular velocity. With no assumption on the halo 
density profile one gets yS M = 0.85, independently of V Q and 
R Q ( |Persic, Salucci 1990| l. However, this is really an absolute 



V(r) = V [l+a (r-R )/R ], 



(8) 



2 We anticipate that p Q is determined here for any reasonable value 
of r o independently of the values taken today for Rd and R e . 
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maximal value and it corresponds, out to R e , to a solid body 
halo profile: V/, oc R" h with 07, = 1. Instead, all mass model- 
ing performed so far for the MW and for external galaxies have 
found a lower value a^Ro) < 0.8, which yields /3 M = 0.77. 
We can also set a lower limit for the disk mass, i.e. p™: first, 
the microlensing optical depth to Baade's Window constrains 
the baryoni c matter within the solar circle to be greater than 
3.9 1O 1O M ( |McMillan, Binney||2009| . Moreover, the MW disk 
B-band luminosity Lb = 2 x 10 lu L o coupled with the very rea- 
sonable value Mr, I Lb = 2 again implies Mo — 4 1O 1O M . AU 
this implies p™ = /3 M /1.3 * 0.65^ We thus take /3 = 0.72!°$ as 
reference range. 

Using the reference values, we get 
p = 0.43^ 1 + 2.9 a Q - 0.64 [fi - 0.72 \ + O.45(r 0D - 3.4 j 



cm J 



-0.l(^-025) + 0.10(,-0.95) 



+ 0.07 



km/ s kpc 



30.3 



(11) 



This equation, which is the main result of our paper, estimates 
the DM density at the Sun's location in an analytic way, in terms 
of the involved observational quantities at their present status 
of knowledge. The equation is written in a form such that, for 
the present reference values of these quantities, the term in the 
square brackets on the r.h.s equals 1, so that the central result is 
p = 0.43 GeV/cm 3 . As such, the determination is ready to ac- 
count for future changes, improved measurement or any choice 
of a B , f3, zo, o), r B D, Q different from the reference values adopted 



here, by simply inserting them in the r.h.s. of eq ( 1 1 



The next step is to estimate the uncertainty in the present de- 
termination of po, which is triggered entirely by the uncertainties 
of the quantities entering the determination. From equation ( fTT) i 
and the allowed range of values discussed above, we see that 
the main sources of uncertainty are a Q ,/3 and r o, which appear 
in the first line. The other parameters give at most variations of 
2-3%, and can be neglected in the following. 

Then, first, it is illustrative to consider a Q , (3 and r QD as inde- 
pendent quantities. We thus have: 



Po = (0. 



43 ± 0.094, 



'(<*o) 



0.016 OT ± 0.096 



\GeV 



CrsD) fcm3 



(12) 



where A( A ) means that A is the total effect due to the possible 
span of the quantity x. 

At this point, we can go one step further, assuming that the 
MW is a typical spiral, and using recent results for the distribu- 
tion of matter in external galaxies, namely that DM halos around 
spirals are self similar (Salucci et al. 20071 and that the frac- 
tional amount of stellar matter f3 shapes the rotation curve slope 
a Q ( |Persic, Salucci|1990| l: 



P = 0.72 - 0.95 or . 



(13) 



3 While these constraints of the disk mass reduce the uncertainty in 
the present determination of p Q , they improve the performance of the 
traditional method very little, where the uncertainties in the disk mass 
value do not trigger the most serious uncertainties of the mass modeling, 
as discussed in the Introduction. 



Using this relation in equation (TT\ we find (neglecting the irrel- 
evant q and zo terms) 



Po = 0.43 



GeV 



cnr 



1 + 3.5 o< + 0.45(r oD -3.4) + 



0.07 



km/s kpc 



-30.3 



(14) 



From the current known uncertainties, with the estimated range 
of a Q , we find 



Po 



0.430 + 0.113 



(ar e ) : 



: 0.096 



(r D) 



\GeV 
/ cm 3 



(15) 



This is our final estimate, which is somewhat higher than pre- 
vious determinations. Its uncertainty mainly reflects our poor 
knowledge of the velocity slope a and the uncertainty in the 
galactocentric Sun distance. 

3. Discussion and conclusion 

In this work we have provided a model-independent kinemati- 
cal determination of p Q . The method proposed here derives p Q 
directly from the solution of the equation of centrifugal equi- 
librium, by estimating the difference between the 'total' density 
and that of the stellar component. 

The method leads to an optimal kinematical determination 
of p , avoiding model-dependent and dubious tasks mandatory 
with the standard method, i.e., a) to assume a particular DM den- 
sity profile and a specific dynamical status for the tracers of the 
gravitational potential, b) to deal with the non-negligible uncer- 
tainties of the global MW kinematics, c) to uniquely disentangle 
the flattish RC into the different bulge/disk/halo components. 

While the measure of p can be performed in an ingenious 
way, it cannot escape the fact that it ultimately depends at least 
on three local quantities, the slope of the circular velocity at the 
Sun, the fraction of its amplitude due to the DM, and the ratio be- 
tween the Sun galactocentric distance and the disk scale-length, 
whose uncertainty unavoidably propagates in the result. 

Two of these three quantities can be related by noting that the 
MW is a typical Spiral and using the relations available for these 
kind of galaxies ( Salucci et al. 2007), so that the final uncertainty 
can be slightly reduced. 

We found that some oblateness of the DM halo and the small 
finite thickness of the stellar disk play a limited role in the mea- 
sure. However, we took them into account by the simple correc- 
tion terms described. 

The resulting local DM density that we find, p = (0.43 ± 
0.11(q. o ) + 0.10( roo )) GeV/cm 3 , is still consistent with previous 
determinations, or slightly higher. However, the determination is 
free from theoretical assumptions and can be easily updated by 
means of equation (Hi as the relevant quantities will become 
better known0 

A final comment is in order. The values of p found in pre- 
vious studies by means of the traditional methods (e.g. Sofue 
|et al. 2009 1 Weber, de Boer|2 009 ) differ among themselves and 
also from the present value only by a small factor. This rela- 
tively good agreement in the values does not imply a concor- 
dance in the underlying mass models, in the various assumptions 
taken or in the data set employed, but is mainly due to the fact 



4 Again, in the traditional method most of the uncertainty in the mea- 
sure of p G discussed in the Introduction cannot be overcome by having 
more and better data. 
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that a> or equivalently A - B (the well known combination of 
Oort constants) is measured with good precision. In fact, from 
Mfi(r) oc V H (f) r we have 



p = *V(1 +2a H ) 



(16) 



where u> 30km/s/kpc, fc 2 is the fraction of the halo contribu- 
tion to the circular velocity at R Q and an = dlogVn/dlogR is 
its unknown slope at the same radius. The quantities k 2 and an 
are experimentally unknown. One can use different assumptions, 
mass modeling and data to get them, but he/she will always find 
that they range from 0.3 (max disk, Persic, Salucci|199 0) to 0.7, 
and from 0.2 to 0.8 (max disk). These are very large variations 
in terms of structural properties of spirals, but only mild ones 
in the determination of p Q , which is dominated by the term u> 2 : 
for any reasonable value of an and k, the density will be in the 
range 0.25 GeV< p < 0.70 GeV. An increase in precision be- 
yond this scale estimate would require an accurate determination 
of k and of the relative mass modeling, of difficult realization by 
means of traditional methods as discussed in the introduction. 
The present method will be able to determine very precisely p Q 
through equation ( 14 1 if improved/new measures of the relevant 
observational quantities also emerge. 



We thus believe that the value given in equation ( 15 i reflects 



the present state-of-the-art knowledge of p and of its uncer- 
tainty, and may result in being very useful in deriving reliable 
future bounds on the DM cross sections involved in direct and 
indirect DM searches. 
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Appendix A: Effect of DM halo oblateness 

In this Appendix we describe the effects of the possible DM halo 
oblateness on our determination of the local DM density. We first 
observe that, while extremal situations (like e.g. a Dark Disk) are 
at odds with observations, a mild oblateness of the DM halo is 
not excluded (see e.g. O'Br ien et al.|2010| for a recent review). 

The effect of a halo oblateness on our method is described 
by the ratio X q (included in equation (|7ji) between the true local 
density of a halo with oblateness q with the one reconstructed 
via the radial derivative Rr 2 d{RV H )ldR: 



PH,q 



(A.l) 



«o,z=0 



where Vn, q is the circular velocity produced by an oblate halo 
PH,q- (Detailed formulae can be found e.g. in Banerjee et al. 



2009 ) Clearly X q reduces to 1 for a spherical halo q — 1, for 
which the spherical Gauss law can be used. In figure A.l we 
have plotted X q as a function of q. As one can see the effect is 
rather weak, of the order of 5%, given a limited oblateness of 
0.9-1, near the upper end. 

It is worth stressing that even a halo oblateness outside 
this reference range can be straightforwardly taken into account 
by the quantity X q . This is another advantage of the proposed 
method with respect to the traditional one. It is in fact much 
easier, to introduce a constant factor of order one in equation 
|7]), than to consider the effect of the halo oblateness in fit- 
ting a global mass model, where it triggers an additional non- 
uniqueness in the structural parameters. 
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Fig. A.l. Effect of the DM halo oblateness q. 



Appendix B: Effect of finite disk thickness 

For the sake of completeness here we also discuss the effect of 
the finite disk thickness. By neglecting it one would slightly 
overestimate the acceleration produced by a given disk mass. 



This can be noted from figure B.l where we plot a correction 
factor for the 'disk' part appearing on the r.h.s. of equation |7}: 



X~ n — 



dr(r 2 a D , Za ) 



d,-(r 2 a Dfi ) 



Re.z=0 



As one can see, given that the disk thickness is measured fairly 
well, z = 0.240-0.250 kpc puric et al.|2008] ) the effect is very 
weak, of the order of 5% (with an uncertainty < 1%). Since this 
is much less than the uncertainty on the disk mass, we neglected 
it in the text, but it may be readily included in the term on the 
r.h.s. of equation (jTjl. 
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Appendix C: Issues in global approaches 

The standard CO method reproduces by a maximum likelihood 
analysis a heterogeneous set of data (among them terminal and 
dispersion velocities, z-motions etc.) with an even more hetero- 
geneous, large set of free parameters (among them R Q , the disk 
and halo masses, the anisotropy in the motions of tracers of the 
gravitational field etc), after which a number of crucial assump- 
tions are taken. In this way it also obtains p Q , as a byproduct. 
This strategy does not guarantee that the global MW model ob- 
tained is physically meaningful, even at the central best fit value 



6 



P. Salucci et al.: Dark Matter density at the Sun's location 



of the parameters. For instance, the global model produces three 
independent and different estimates of V(r): 1) from the HI ter- 
minal velocities WW (via V(r) = Vr(r) + V Q r/R Q ), 2) from the 
halo stars dispersion velocities (via the Jeans equation), and 3) 
from the gravitational potential produced by the MW mass dis- 
tribution. Clearly, by physical consistency, these three estimates 
of V(r) must agree, but in the strategy commonly employed there 
is nothing forcing this to happen. 

In fact, let us look at the galaxy global model obtained re- 
cently by this method by Catena, Ullio (2009). The above 3 dif- 



ferent estimates of V(r) are obtained by using the best fit model 



parameters as given in their Tables 2 and 3. We find (figure C.l 



that they disagree by 10% in amplitude and 0.3 in slope, quanti- 
ties more than allowed by the observational errors. This implies 
that this method has an intrinsic uncertainty that may lead to a 
biased measure of p Q and of its uncertainty. 
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Fig. C.l. Different estimates of the circular velocity of the Milky 
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and velocity dispersions ( Xue et al.|2008] > (black points). 
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